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Abstract 

In this paper we analyze the vacuum expectation values of the field squared and the 
energy-momentum tensor associated to a massive scalar field in a higher dimensional cosmic 
string spacetime, obeying Dirichlet or Neumann boundary conditions on the surface orthog- 
onal to the string. In order to develop this analysis the corresponding Green function is 
obtained. The Green function is given by the sum of two expressions: the first one corre- 
sponds to the standard Green function in the boundary-free cosmic string spacetime and the 
second contribution is induced by the boundary. The boundary induced parts have oppo- 
site signs for Dirichlet and Neumann scalars. Because the analysis of vacuum polarization 
effects in the boundary-free cosmic string spacetime have been developed in the literature, 
here we are mainly interested in the calculations of the effects induced by the boundary. In 
this way closed expressions for the corresponding expectation values are provided, as well 
as their asymptotic behavior in different limiting regions is investigated. We show that the 
non-trivial topology due to the cosmic string enhances the boundary induced vacuum polar- 
ization effects for both field squared and the energy-momentum tensor, compared to the case 
of a boundary in Minkowski spacetime. The presence of the cosmic string induces non-zero 
stress along the direction normal to the boundary. The corresponding vacuum force acting 
on the boundary is investigated. 



PACS numbers: 98.80. Cq, ll.10.Gh, 11.27.+d 



1 Introduction 

Cosmic strings are topologically stable gravitational defects which may have been created in 
the early Universe after Planck time by a vacuum phase transition [H [2]. The gravitational 
field produced by a cosmic string may be approximated by a planar angle deficit in the two- 
dimensional sub-space. The simplest theoretical model which describes a straight and infinitely 
long cosmic string is given by a Dirac-delta type distribution for the energy-momentum tensor 
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along the linear defect. Also this object can be described by classical field theory, coupling 
the energy-momentum tensor associated with the Maxwell-Higgs system investigated by Nielsen 
and Olesen in [3] with the Einstein equations [U [5]. Though the recent observational data 
on the cosmic microwave background have ruled out cosmic strings as the primary source for 
primordial density perturbation, they are still candidate for the generation of a number of 
interesting physical effects such as gamma ray bursts [6] , gravitational waves [7] and high energy 
cosmic rays [8]. Moreover, in the framework of brane inflation [9]-[TT], cosmic strings have 
attracted renewed interest, partly because a variant of their formation mechanism is proposed. 

Although the geometry of the spacetime produced by an idealized cosmic string is locally flat, 
the planar angle deficit provides nonzero vacuum expectation values (VEVs) for different physical 
observables. In this context, the VEVs of the energy-momentum tensor have been calculated 
for scalar and fermionic fields in |12j-[16j and |17]-[20j. respectively. Another type of vacuum 
polarization takes place when boundaries are present. In this sense, by imposing boundary 
conditions on quantum fields, additional shifts in the VEVs of physical quantities, such as the 
energy density and stresses, take place. This is the well-known Casimir effect (for a review 
see [21]). The analysis of the Casimir effect in the idealized cosmic string spacetime have been 
developed for scalar [221. vector [231 IS] and fermionic fields [25j[26] , obeying boundary conditions 
on cylindrical surfaceso The Casimir force for massless scalar fields subject to Dirichlet and 
Neumann boundary conditions in the setting of the conical piston has been recently investigated 
in |28| . Continuing along this line of investigation, in this paper we shall analyze the contribution 
on the vacuum polarization effects in a higher-dimensional cosmic string spacetime induced by 
a scalar field obeying Dirichlet or Neumann boundary conditions on a surface orthogonal to 
the string. In addition to be a new perspective related with the Casimir effect, this present 
investigation may be relevant in the analysis of vacuum polarization effects induced by a brane 
in anti-de Sitter spacetime. 

The paper is organized as follows. In section [2] we provide the general expression for the 
scalar Green function in a higher-dimensional cosmic string spacetime admitting that the field 
obeys Dirichlet or Neumann conditions on a boundary orthogonal to the string. We shall see 
that this Green function is expressed in terms of two distinct contributions. The first one is the 
standard Green function for a massive scalar field in a boundary-free cosmic string spacetime, 
and the second contribution is due to the boundary condition obeyed by the field operator. The 
first contribution to the Green function is divergent at the coincidence limit; as to the second 
one, it is finite in this limit for points away from the boundary. Moreover, for specific values 
of the parameter associated with the planar angle deficit, the complete Green function can be 
expressed in a closed form in terms of a finite sum of the Macdonald functions. Because the 
analysis of the VEVs of physical quantities induced by the cosmic string have been developed 
in literature by many authors, in sections [3] and E] we shall calculate the contributions to the 
VEVs of the field squared and the energy-momentum tensor induced by the boundary. We shall 
see that near the boundary and for points outside the string these contributions become more 
relevant than the contributions due to the cosmic string itself. In section [5j we summarize the 
most important results obtained. In this paper we shall use the units H = c = 1. 

2 Green function 

2.1 Geometry of the problem and the heat kernel 

We consider a massive scalar quantum field propagating in a D-dimensional cosmic string space- 
time. By using the generalized cylindrical coordinates with the cosmic string on the subspace 

1 Also vacuum polarization effects induced by a composite topological defect have been analyzed in |27| . 
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defined by r = 0, being r ^ the radial polar coordinate, the corresponding metric tensor is 
defined by the line element below: 

D-l 

ds 2 = g ik dx i dx k = -dt 2 + dr 2 + a 2 r 2 d(p 2 + dz 2 + V] (dx') 2 . (1) 

The coordinate system reads: x % = (t,r,ip, z,x l ), with ip E [0, 27r], and t, z, x l G (— oo, oo). 
The parameter a, smaller than unity, codifies the presence of the string. In a 4-dimensional 
spacetime, this parameter is related to the linear mass density of the string by a = 1 — 4G/x, 
with G being the Newton gravitational constant. In this analysis we shall admit the presence of 
extra coordinates, x , defined in a Euclidean (D — 4)-dimensional subspace. 

For a scalar field propagating in an arbitrary curved spacetime the field equation has the 
form: 

(□ - m 2 - £R) 4>[x) = , (2) 

with □ denoting the covariant d'Alembertian and R is the scalar curvature. In ([2]) we have 
introduced an arbitrary curvature coupling £. The minimal coupling corresponds to £ = and 
for the conformal one £ = £ c = (D — 2)/[4(D — 1)]. We shall assume that the field obeys Dirichlet 
boundary condition on the hypersurface orthogonal to the string and located at z = 0: 

(j>(x) = o,2 = 0. (3) 

The Green function associated with a massive scalar field in a curved spacetime obeys the 
second order differential equation 

(□ _ m 2 _ £ R j G ^ x / ) = _ S D^ x>) = J D {X-X') ^ (4) 

where S D (x, x') represents the bidensity Dirac distribution. This Green function can be obtained 
within the framework of the Schwinger-DeWitt formalism as follows: 

POO 

G(x,x') = / ds K(x,x'; s) , (5) 
J o 

where the heat kernel, lC(x,x';s), can be expressed in terms of a complete set of normalized 
eigenfunctions of the operator defined in as follows: 

K{x, x'; s) = J2 $*(x)K(x')e- Sa2 , (6) 

17 

with a 2 being the corresponding positively defined eigenvalue. 
Writing 

[n-m 2 -tR)$ c {x) = -a 2 $ a {x) , (7) 

in the spacetime defined by the line element ([T|), a complete set of normalized solutions of the 
above equation, compatible with the boundary condition ([3]), can be specified in terms of a set 
of quantum numbers (w, q, n, k z , fc/), where n = 0, dbl, ±2, ... , (oj, k[) G (— oo, oo) and 
(<?; k z ) ^ 0. These functions are given by: 

e i(n^>+k-x— wt) 

* ff (x) = 2 y~ ( 27r )(g-i)/2 J i"i/^(g r ) sin ( A: ^) ' ( 8 ) 
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being J v (x) the Bessel function and x = (x 4 , . . . ,x D 1 ). The corresponding positively defined 
eigenvalue is given below: 

a 2 = J 2 + q 2 + k 2 z + k 2 + m 2 . (9) 

The heat kernel can be given in terms of the above eigenfunctions according to ([6]). After 
performing the integrals with the help of [29], we obtain: 

2p~^~ sm2 ( zz'\ a frr'\ 

K{X ^ S) = l<4^ Sinh (27) £ ' "'HA" 27 ' <10) 

where is the modified Bessel function and 

Ap 2 = -At 2 + r 2 + r' 2 + Ax 2 + z 2 + z' 2 , (11) 

with A</? = ip — tp' , At = t — t', Ax = x — x'. 

In general, it is not possible to provide a closed expression for the Green function by inte- 
grating over the variable s the heat kernel function (jlOp . according to ©. However, for massless 
fields and for specific values of the parameter a, the corresponding Green functions can be 
expressed in terms of a finite sum of the associated Legendre functions and the Macdonald 
ones, respectively. These two different situations will be analyzed separately in the following 
subsections. 

The case of a scalar field with Neumann boundary condition, d z (j) = at z = 0, can be 
considered in a similar way. The corresponding eigenfunctions have the form ([8]) with the 
replacement •svrt(k z z) —> cos(k z z). The expression for the heat kernel is obtained from ()10p with 
the replacement sinh(zz'/(2s)) — > cosh(zz' / (2s)) . 



2.2 Special case 

The analysis of vacuum polarization effects associated with a quantum scalar field in a cosmic 
string spacetime have been developed by many authors for the case where the parameter a is 
equal to the inverse of an integer number p, i.e., when a = 1/p (see p3]-[15]). In this case the 
corresponding Green function can be expressed in terms of p images of the Minkowski spacetime 
function. Recently the image method was also used in [22] to obtain a closed expressions for 
massive scalar Green functions in a higher-dimensional cosmic string space-time obeying Robin 
boundary condition on a cylindrical surface coaxial with the string. Here, in this subsection, we 
shall consider this specific situation, i.e., a being the inverse of an integer number, to obtain the 
Green function in a closed form for the physical situation under consideration. For this case, 
the expression of the heat kernel can be further simplified with the help of the formula [30|, 13 lj : 

+00 ^ p— 1 , 

£ e^I ln{p (rr'/2s) = ± • (12) 

n=—oo P k=0 

The corresponding heat kernel reads, 

/C(x,x^) = ^^sinh(^)ge-^ , (13) 

where 

V fc = -At 2 + Ax 2 + z 2 + z' 2 + r 2 + r' 2 - 2rr' cos (Aip/p + 2-rrk/p) . (14) 
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Finally, substituting the above function into with the help of we get, 

G ( x ' x ') = ?7 \D/2 Yl [/-D/2-1 ( mV H-)) ~ fo/2-i (mV k{+) )] , (15) 

where 

V fc(=F) = [-At 2 + Ax 2 + (z =F .z') 2 + r 2 + r' 2 - 2rr' cos (A99/P + 2irk/p)} 1/2 . (16) 

In (I15p and in what follows we use the notation 

/„(x) = K v (x)/af, (17) 

being K u (x) the Macdonald function. The expression (|15p is further simplified in the case of 
a massless field. By using the asymptotic of the Macdonald function for small values of the 
argument [32], one finds 

r(D/2-l)^r ._ D „ 



47r D/2 

fc=0 

We can see that the Green function (|15p vanishes for 2 or z' being equal to zero. It can be 
presented as the sum of two different contributions as shown below: 

G(x, x') = G cs {x, x') + G h (x, x') . (19) 

The first term in the right-hand side coincides with the Green function for a massive scalar 
field in the absence of the boundary. It is divergent at the coincidence limit and the divergence 
comes from the k = term. As to the second contribution, Gb(x,x'), it is a consequence of the 
boundary condition imposed on the field. This contribution is finite at the coincidence limit for 
points outside the boundary. 

The renormalized Green function, used for the evaluation of finite and well defined VEVs, 
is given by subtracting from the corresponding Green function the Hadamard one. Because the 
cosmic string spacetime is locally flat, the Hadamard function coincides with the Green function 
in the Minkowski spacetime. For the above Green function the renormalization procedure can 
be implemented explicitly by discarding the k = component of the function G cs (x,x'). 

The formulae for the Green function in the case of Neumann boundary condition are obtained 
from (fl"5|) and (fTSj) changing the sign of the terms with V&( + ). Consequently, the boundary 
induced parts in the Green function, G\>(x,x'), for Dirichlet and Neumann scalars differ by the 
sign. This is the case also for the corresponding parts in the VEVs of the field squared and the 
energy- momentum tensor discussed below in section [3] and [H 

2.3 General case 

For general case where p = 1/a is not an integer number, the Green function can be expressed in 
an integral form by substituting the heat kernel (jlOl) into ([5]) . After some intermediate steps the 
Green function is presented in the form (|19p . where the boundary- free and boundary induced 
parts are given by the expressions 

G cs (x,x') = £ e m M dw wW-te-P"-^ J |n|p (rr'w/2) , 

G ^ X ') = -(4^572 E ^jj, dww^e-^-t I Hp (rr'w/2) , (20) 
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where 



(T) 



-At 2 + Ax 2 + (z =F ) + r 2 + r 



,,'2 



(21) 



We can also verify that the Green function vanishes for z or z' equal to zero. Here we can analyze 



the behavior of the both contributions in the coincidence limit. In this limit, V, 



(-) 



2r 2 so 



due to the exponential behavior of the modified Bessel function for large arguments, the first 
integral diverges for large values of the variable w; as to the second contribution, V( + ) = 4z 2 +2r 2 , 
consequently for z / the integrand in the expression for Gb(x,x') goes to zero exponentially 
for large values of the integration variable. 

We can also provide a closed expression for the Green function in the limit of a massless 
field. With the help of [29] we find: 



G(x,x') 



-i(D-3)7r/2 



P 



(2vr)( 13 + 1 )/ 2 (rr 1 ) / 2 - 1 



£ ■ 

l= — oo 



Q\n\p-l/l( COshu (-)) 



Q\n\ P -l/l( COshu (+h 



(sinh u H ) ( D - 3 )/ 2 (sinh u {+) ) ( D ~ 3 )/ 2 



(22) 



being Qu(x) the associated Legendre function and 



cosh-u 



(=F) 



(=F) 



2rr' 



> 1 



(23) 



Finally, we have to say that the renormalized Green function can be obtained by subtracting 
from the Green function the corresponding function in the Minkowski spacetime, which is given 
by G cs (x, x') taking p = 1. 



3 Vacuum expectation value of the field squared 

This and the following sections will be devoted to the calculations of vacuum polarizations 
effects induced by the boundary. Two main calculations will be performed. The evaluation of 
the VEV of the field squared, in the first place, followed by the evaluation of the VEV of the 
energy- momentum tensor. 

The VEV of the field squared is formally given by evaluating the Green function at the 
coincidence limit. In this analysis the complete Green function is given by the sum of the 
Green function in the cosmic string spacetime in the absence of the boundary plus the boundary 
induced part. In this way we may write, 

<0 2 > = <</> 2 > cs + (0 2 ) b . (24) 

However, because the singular behavior of G cs (x, x') at the coincidence limit, the renormalization 
procedure is needed for the first contribution of the above expression. The second contribution 
is finite at the coincidence limit for points outside the hypersurface z = 0. Because the VEV of 
the field squared in the cosmic string spacetime has been analyzed by many authors, here we 
are mainly interested in the analysis of the quantum effects induced by the boundary. 

According to the previous section, we shall analyze the VEV of the field squared induced 
by the boundary for p being an integer number in the first part, and for the general case in the 
second one. 
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3.1 Special case 

Being p an integer number, the VEV of the field squared induced by the boundary is given 
simply by taking the coincidence limit of G\>(x',x) given in f)15|) . The result is presented by the 
expression 

,fl-2 P- 1 



s.2\ 



m 



>b 



&|(V2 E/^(24 2 +A|), (25) 
where 

Sk = sm(irk/p) . (26) 

As it is seen, the boundary induced part in the VEV is negative. In (|25|) . the part with the term 
k = is the corresponding VEV in the Minkowski spacetime and we can write: 

(0 2 >b = <^> b P=1) + <<^> b C , (27) 

where the second term on the right-hand side is the part of the VEV induced by the non-trivial 
topology of the cosmic string spacetime. Note that for the renormalized pure topological part 
one has the expression (see [22 ) 

m D - 2 ^ „ 

(28) 



III . 

-^7572 2^ fo/2-i (2mrs k ) . 



The latter is always positive. It is of interest to note that the topological part in the total VEV, 
(</ )2 )cs + (^ 2 )b ) vanishes on the boundary. 

For z / and for points far from the string, r 3> \z\, the dominant contribution to the 
boundary induced part (j25|) comes from the k = term and to the leading order we get 



(<P\ « (^ 2 )[f- 1) = (2m\z\) . (29) 

Note that for z ^ and for points on the string axis, r = 0, one has (</> 2 )b,r=o = Pi^ 2 )^ -1 ^ ■ I n 
the limit \z\ ^> r and m\z\ 3> 1, the leading order term provides an exponentially suppressed 
behavior below, 

2 p m (D-3)/2 e -2m\z\ 

For a massless field, from (|25p we obtain 

Ia2\ F(D/2 — 1) ^ 2 . 2 2\l-D/2 /qi \ 

^) b = ( 47r) D/2 1J Z + r S k) ■ ( 31 ) 

In this case, at large distances from the boundary, \z\ 3> r, the boundary induced VEV decays 
as power-law: (0 2 )b °c r 2 ~ D . In figure [1] we exhibit the behavior of (|25p as a function of the 
dimensionless variables mr and mz for D = 4, p = 3. 

As we have already explained and it is seen from the graph, for fixed z the boundary induced 
contribution is finite for z ^ and r = 0. On the other hand, for fixed non-vanishing radial 
coordinate, near the boundary it is dominated by the k = term (see (|29|) ) . When z goes to 
infinity the VEV is exponentially suppressed. 

In figure [2] we exhibit the behavior of (|25h in the case of a D = 4 cosmic string as a function 
of the dimensionless variable mr for three distinct values of p = 2, 3, 4 and for mz = 0.5. We 
can see that the effects induced by the string become more relevant for larger values of p. 
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Figure 1: The boundary induced part in the VEV of the field squared, (<f, 
dimensional cosmic string spacetime with function of mr and mz. 
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Figure 2: This graph provides the behavior of (</> 2 )b/?7i 2 in a 4-dimensional spacetime as a 
function of mr for mz = 0.5 and for several values of the parameter p (numbers near the 
curves) . 
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3.2 General case 



In the case of p being not an integer number, the VEV of the field squared induced by the 
boundary is obtained by taking the coincidence limit of Gb(x',x) given in (|20p . After a change 
of the integration variable it can be written as: 

2—D roc 

(<t>\ = -^yJJ- 2 l dy ^-2^(2*7^1)^7(2,) Sp{y) (32) 

with 

oo 

s p(y) = T np(y) ' ( 33 ) 

n=0 

where the prime means that the term with n = should be halved. As it is seen from (|32p . 
the boundary induced VEV is negative for general case of p. Although being not possible to 
provide a closed expression for the above VEV, some limiting cases can be obtained. Due to the 
exponential decay, for \z\ ^> r the dominant contribution comes from the region near the lower 
limit of the integration. In this region we may approximate S p (y) ~ 1/2 and with the help of 

[29] we obtain {4> 2 }b ~ p(<^ 2 )b ^ • ^ e can see ^ na * f° r P > 1 the presence of the cosmic string 
increases the above results when compared with the corresponding ones in the absence of it. 
We can find an integral representation for the VEV of the field squared by using the formula 

S M = I £' - *!M r ix (34) 

P ^ 7r Jo cosn(px) — cos(p7r) 

where po is an integer number defined by 2po < p < 2pQ + 2. For even values of p the corre- 
sponding formula is obtained from (|34l) by taking the limit. In this limit the second term on the 
right-hand side becomes e~ y /(2p). The formula (|34p is obtained as a special case of the more 
general formula derived in |26j . 

Substituting (|34p into (|32p , the integration over y is done explicitly in terms of the Macdonald 
function and one finds 



2m D - 2 



Po 



(2tt) d / 2 

P ■ i „\ a fD/2-i(^ m V z 2 + r 2 cosh 
cosh(2px) — cos(p7r) 



-—sm(jm) / dx 
v Jo 



k=0 

f n m i (2m,\/ z 2 -I- r 2 cosh 2 x) 

(35) 



It can be seen that for integer values p this result reduces to the expression (J25J). On the string 
axis for the boundary induced VEV one has 

(<p\,r=o=p(^) { r 1) - (36) 

For points near the boundary, \z\ m~ 1 , \z\ <C r, the main contribution to the boundary 
induced VEV comes from the term k = and to the leading order we find 

2 ^ rp/2-1) 

{4> >b ~ -^)D/ M D-2- ( 37 ) 

This leading term does not depend on the angle deficit and it coincides with the corresponding 
term for a boundary in Minkowski spacetime. At large distances from the boundary, \z\ 3> m _1 , 
the boundary induced part is exponentially suppressed and the VEV is dominated by the pure 
topological part (4> 2 ) cs . 
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For a massless field, by using the relation f u (y) ~ 2 V T(v)y 



,-2v 



0, we find 



2T(D/2 - 1) 
(4vr)^/2 



Po 
fc=0 



^ + r 2 s 2)l-D/2 



■— sin(p7r) 

7T 



1-13/2 -. 



00 (z 2 + r 2 cosh^ 
o cosh(2px) — cos(p7r) 



(38) 



In this case, at large distances from the boundary, \z\ 3> r, the leading term in the corresponding 
asymptotic expansion has the form 

P T{D/2 - 1) 



■ v(6 2 ) {p=1) - 



(39) 



Recall that near the boundary we have the behavior given by (|37p . 

For a scalar field with Neumann boundary condition at z = 0, the corresponding formulae 
for the VEV of the field squared are obtained from those given above by changing the sign of 
the boundary induced part. 

In figure [3] we plot ((p 2 )^/m 2 as a function of the parameter p for mz = 0.5 and for several 
values of mr (numbers near the curves) in a 4-dimensional spacetime (D = 4). Note that for 
the first term in the right-hand side of (I27D . for mz = 0.5 one has {4> 2 )u ~ —0.0152m 2 . 
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Figure 3: The topological part in the boundary induced VEV of the field squared, ((/> 2 )^/m 2 , in 
a 4-dimensional cosmic string spacetime as a function of p for mz = 0.5 and for several values 
of mr (numbers near the curves). 

We can also analyze the topological part in the boundary induced VEV of the field squared 
by adding and subtracting from (|32j) the part corresponding to the Minkowski spacetime (p = 1): 



„2-D 



(2n)D/2 J 
J2-D 



(2tt) d / 2 



X3 

dyy N/2 e -(2zyr^l)y- m ^/(2y) ^ ^y) 

n=0 
oo 

d yym e -^y + l)y-m^/{2y) £ ^Wtf) _ . (4Q) 



n=0 



E/oo 
I n (y) = e y 12. 
71=0 

Of course, the corresponding expression coincides with (|29p which is independent of r and 
diverges for z = 0. 
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An alternative representation for the topological part can be obtained by using the Abel- 
Plana formula (see, for instance, [33]) for the summation over n in the second term on the 
right-hand side of ([3D]): 

v^'w \ f°° , „, x f 00 , F(iu)-F(-iu) 

^F(n) = duF(u) + i du K > _\ ' . (41) 



ra=0 



In our case, F(n) = pl np (y). Now we can see that in the evaluation of the difference the terms 
coming from the first integral on the right-hand side of the Abel-Plana formula cancel out and 
one obtains 

00 / 2 r°° 

Y\ Wn P (y) - I n (y)) = ~ dv g(v,p)K iv {y) , (42) 
^0 nJ ° 

where we have introduced the notation 

g(v,p) = sinh(^) ( — j- - } ) . (43) 



e 2irv/p _l g 



The respective contribution becomes: 

D foe 

dyy 



. 2 D poo 

2\C = 4r / -(2 2 2 /r 2 +l) J/ -m 2 r 2 /(2 2 /) 



b (2^)^/2+1 J Q 

/>oo 

x dv g{v,p)K iv (y) . (44) 
Jo 

For r/0 the above correction is finite at the boundary. 

Taking z = 0, for a massless field the integral over y in (|44[) is evaluated in terms of the 
gamma function and we obtain 

a(A 7r \-{D+l)/2 poo 

(^Xz=o = - T{{D -l)/2)rO-2 J, dV 9{V ' P) |r {D/2 ~ 1 " W)? • (45) 

For an even dimensional spacetime the modulus of the gamma function in this formula is ex- 
pressed in terms of elementary functions and the integral can be explicitly evaluated. In par- 
ticular, for D = 4 one has |T(1 — iv)\ 2 = wv/ sinh(7ro) and the integral above becomes equal to 
tt(p 2 — l)/24. In this case we get: 



,2 )C _ P -1 



(46) 



This result is also obtained from (|25p . in the special case of integer p. Taking D = 4, for a 
massless fields we have: 



p-i 



l07T z z z + H sin (ttk/p) 

On the boundary, z = 0, the sum in this formula is explicitly evaluated, Y^k=i s ^ n _2 ( 7r ^/p) = 
(p 2 - l)/3, and (g7} reduces to (gH). 
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4 Energy- momentum tensor 

Following the same line of investigation, in this section we are interested to calculate the con- 
tribution induced by the boundary on the VEV of the energy- momentum tensor. Similar to the 
case of the field squared, the energy-momentum tensor is presented in the decomposed form, 



(Tik) — (Tik) cs + (Tik)bi 



(48) 



where {Ti k ) cs corresponds to the geometry of the cosmic string without boundaries. In order to 
evaluate the boundary induced part we shall use the following expression: 



(T ik ) h = lim di,d k G b (x, x') + [(£ - 1/4) g ik D - £V,V fc - £R ik ] (0 2 ) b , 



(49) 



where for the spacetime under consideration the Ricci tensor, Ri k , vanishes. 

By using the expression f)35f) . for the covariant d'Alembertian appearing in (|49|) one finds 



□(0 2 )b 



8m 



D 



{2tt) d / 2 

p sin (jot) 
2~k 



VQ 



where 



dx 
it'i 

Wo 



l±m\z 2 + str 2 )f D/2+1 ( Wl ) - (1 + 2s 2 )f D/2 ( Wl )] 

1=0 

w\ + m 2 r 2 sinh 2 x) f D / 2 +i (w 2 ) - (2 + cosh x)f D / 2 (w 2 ) 



cosh(px) — cos(p7r) 



. (50) 



2m \ z 2 + r 2 si, 



2mJ z 2 + r 2 cosh (x/2) 



(51) 



In the evaluation of the derivative term d^d^G^x^x') we need the expression for the sum 
^2n=-oo n2 I\n\p (y)- For the summation of this series we use the relation 



J2 n 2 I ln[p (y) = 2p~ 2 (y 2 d 2 + yd y - y 2 ) S p (y), 



(52) 



and the formula (|34p . In the evaluation of the other derivative terms we can put <j> = (f)' before 
the differentiation. In this case, by using the formula (|34p . for the boundary induced part in the 
Green function one finds 



G b (x,x') 



2m 



D-2 



llp'=lp 



(2ir) D / 2 
psin(p7r) 

2~K 



fD/2~i{m\JV{+) ~ 2rr'cos(27rl/p)) 
fD/2-i(my/V(+) + 2rr'cosha;' 



1=0 

oo 



dx- 



cosh(px) — cos(p-7r) 



(53) 



where V( + ) is defined by ([211 . 

Combining the formulae given above, for the separate components of the energy-momentum 
tensor we find the following VEVs: 



2m 



D 



(2tt) d / 2 



Vt) 



p 

7T 



s'm(pir) 



F-(2mr,2mz,s k ) 

k=0 

00 F- (2mr, 2mz, cosh x) 
o cosh(2px) — cos(p7r) 



(54) 
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Here the functions are defined as 

F°(u,v,w) = ^-l)(w 4 u 2 + v 2 )f D/2+1 (y)+[l-(4C-l)(l + 2w 2 )]f D/2 (y), 

F}(u,v,w) = (^-l)v 2 f D/2+l (y)+2[l-2ai + w 2 )]f D/2 (y), 

F 2 (u,v,w) = [4^w i u 2 + v 2 )-y 2 ]f D/2+1 (y) + 2[l-2^1 + w 2 )}f D/2 (y), (55) 

Fl(u,v,w) = (4£ - l)w 2 [w 2 u 2 f D/2+1 (y) - 2f D/2 (y)] , 

Fi(u,v,w) = -(4£ - l)uvw 2 f D/2+1 (y), 

the indices 0, 1, 2, 3 correspond to the coordinates t, r, ip, z, and 



y = \J v 2 + u 2 w 2 . (56) 

For the components (no summation over I) (T/)b, I = 4, . . . , D — 1, one has the relation (T/)b = 
(Tg)b- Of course, the latter relation is a direct consequence of the invariance of the problem 
with respect to the boosts along the directions x , I = 4, . . . , D — 1. Due to the presence of the 
boundary the boost invariance along the z-axis is lost and, as a consequence, (Tf )b / {Tq)\>. 
As we see, the vacuum energy-momentum tensor is non-diagonal. Note that the off-diagonal 
component vanishes for z = and r = 0, separately. In addition, on the string axis, r = 0, 
one has (T 2 2 )b = (^i)b- Eor a scalar field with Neumann boundary condition at z = 0, the 
corresponding formula for the VEV of the energy-momentum tensor is obtained from (I54D by 
changing the sign of the boundary induced part. 

Due to the presence of the off-diagonal component, from the covariant conservation condition, 
V A j(T ! ^) = 0, for the energy-momentum tensor two non-trivial differential equations follow: 

d r (r(T?) h ) + rd z (T?) h = (T 2 ) h (57) 

and 

0,<2f>b = -ld r (r(Ti) h ) . (58) 

It can be checked that these relations are obeyed by the VEV of the energy-momentum tensor 
above. 

The k = term in (|54p is the corresponding VEV for a flat boundary in Minkowski spacetime. 
The corresponding expression takes the form (no summation over I) 

(Tbt^ = -J0^D]2 ^ ~ l) (^z) 2 f D/2+l (2m\z\) + 2 (1 - 20 f D/2 (2m\z\)] , (59) 

for I = 0, 1, 2, 4, . . . , D — 1, and (T^)^ - = for the other components. Note that for a 

conformally coupled massless field (^) b P_1 ' ) = 0. 

For points near the boundary, \z\ <C m~ l , \z\ <C r, one has 

(Ti) b « (Tf)^ - ~2(D - 1) (e - ^P^, (60) 

for / = 0, 1, 2,4, . . . , D — 1, and these components diverge on the boundary. For a conformally 
coupled field the leading term vanishes and the subleading term should be taken: 

^Wb-^Wb - (47r)0/2 (£)_!)(£) - 2 )' lD j 

The corresponding energy density, (T 00 )b = — (Tg )b, is negative for both minimally and confor- 
mally coupled scalar fields. For the components (Tf )b, (T^ )b> the = term in (|54p vanishes 
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and these components are finite for r / 0. Note that for z / the boundary induced part in 
the VEVs is finite on the string axis, r = 0. By taking into account that the pure topological 
part diverges on the string we conclude that the latter dominates for points near the string. 

At large distances from the string, r 3> m^ 1 , r S> \z\, the dominant contribution to the 
boundary induced VEV (|54p comes from the k = term and, to the leading order, the VEVs 
coincide with the corresponding expressions for a boundary in Minkowski spacetime. In this limit 
the effects due to the non-trivial topology of the cosmic string are exponentially suppressed. At 
large distances from the boundary, \z\ 3> m _1 , the boundary induced VEVs in the components 
of the energy-momentum tensor are suppressed by the factor e _2m l 2 l. 

In the case of a massless field, for the VEVs we find the expression 



P 

7T 



rp/2) 

(4vr) D / 2 



sin(p7r) 



£' 

fc=0 



dx 



F\ 0)i {r,z,s k ) 
(z 2 + r^l) /^ 1 
F[ Q )i(r,z, cosh x) 



cosh(2px) — cos(j>7r) (z 2 + r 2 cosh 2 x) D / 2+1 



(62) 



The functions for the separate components are defined by the expressions 



(0)0 



[r,z,w) 



{(4£ - 1)[(D - 2)w 2 - 1] + 1} r 2 w 2 + [(4£ - 1)(D - 1 



2w 2 ) + ll z 2 , 



F} 0)1 (r,z,w) = [2-4£(l + w 2 )] r 2 w 2 + [4£ (D - 1 - w 2 ) - D + 2} z 2 , 



Ff Q)3 (r,z,w) 



[4£ {{D - l)w 2 - 1) - D + 2] r 2 w 2 + [4£ (D - 1 
(4f - l)w 2 [{D - 2)w 2 r 2 - 2z 2 ] , 
-D(4£ - l)rzw 2 . 



w 



D + 2]z 2 , (63) 



Now it can be checked that for a massless conformally coupled scalar field the tensor (T-)^ is 

traceless. Note that for a massless field the VEV (^i)jf is given by the right-hand side of 
(j60|) and vanishes for a conformally coupled field. The expression on the right-hand side of (j60|) 
gives the leading term in the asymptotic expansion of (Tf)b at large distances from the string, 
r ^> \z\, for the components with i = I = 0, 1, 2, 4, ... ,D — 1. In the same limit, the other 



non-zero components decay as {T^)\, oc r D and {T^)^ oc r 



-D-l 



At large distances from the 



boundary, \z\ 3> r, the boundary induced part in the VEVs of the diagonal components behave 
as (no summation over I) (T/)b °c whereas for the off-diagonal component one has (T^)b 



oc z 



-D-l 



In the case p being an integer number, the general formula (J54J) takes the form 

p-i 



m 



D 



(2tt) d / 2 



} F-(2mr, 2mz,s k ), 



k=0 



(64) 



with the functions F-(u,v,w) defined in (|55p . Note that the corresponding expression for the 
boundary-free part has the form [22] (no summation over I) 



07) 



in 



D 



p-i 
k=l 



where 



F csfi (u,w) 
F CBil (u,w) 



{2ir) D / 

(1 - 40u 2 w 4 f D/2+l (uw) - [1 + 2(1 - 40- 2 
(4£ 



(65) 



/d/2( 



UW 



w 



(4e 



l) fD/2{uw), 
!) [/d/2(^) - « 2 ^ 2 /d/2+i(^)] , 



(66) 
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and Fes K u > w ) = Fcs,o( u i w ) f° r / = 3, . . . , JD — 1. For a massless field the formulae 
are simplified to (no summation over I) 



and 



T{D/2)5[ P ^ F^(s k ) 



2 ( 47r )D/2 r D 



k=l 



£CP/2) 

2(4vr)^/ 2 ^ (z 2 + r 2 4) D / 2 + 1 ' 



E 



(67) 
(68) 



with the functions 



(D - 2) (1 - 4£)w 2 - 1, I = 0, 3, . . . , D - 1, 



1 - D 



A^w 2 



(69) 



and the functions F^(r, z, w) are given by (|63p . For a conformally coupled massless scalar field 
the vacuum energy density (T 00 ) is always positive. 

In the left / right panel of figure H] the VEV of the energy density induced by the boundary is 
exhibited for a minimally /conformally coupled massive scalar field in a 4-dimensional spacetime, 
as a function of mr and mz for p = 3. We can see that the energy density crucially depends on 
the curvature coupling parameter £. 





Figure 4: The boundary induced part in the VEV of the energy density, (T 00 )b/m 4 , for minimally 
(left plot) and conformally (right plot) coupled scalar fields in a 4-dimensional spacetime as a 
function of mr and mz for p = 3. 

The normal vacuum force acting on the boundary is finite for r / and is determined by 
the component (T| )b evaluated at z = 0: 



(Ti)b,z=o 



2m u (l - 4f) 
(2tt) d / 2 



pa 



^2s 2 k F 3 (2mrs k ) 



k=l 



-— sin(p7r) 

7T 



00 F 3 (2mr cosh x) cosh 2 x 
cosh(2px) — cos(p7r) 



where we have defined the function 



F 3 (u) = u 2 f D/2+1 (u) - 2f D/2 (u). 



(70) 
(71) 
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The vacuum effective pressure on the boundary is given by P = (T 3 3 ) b)2=0 . Note that the 
dependence of the force on the curvature coupling parameter appears in the form of the factor 
(1 —44;). For a massless field, the expression ([7U|) reduces to 



(?f)b,z=0 



(D 



(1 - 4£)r( J D/2) 



(4 7r )D/2 r D 



■— sin(p7r) 

7T 



dx- 



Po 

fc=i 

cosh 2 " 13 



D 



X 



i cosh(2p:r) — cos(p7r) 
In particular, for p being an integer number and for D = 4 one has 

1-44 



(T 3 )b,*=o 



48vr 2 r 4 



(P 2 



1) 



(72) 



(73) 



For the both minimally and conformally coupled scalar fields the corresponding effective pressure 
is positive. In the left plot of figure [5] we have presented (T| )b,z=o/'m D for a minimally coupled 
scalar field in D = 4 as a function of mr for separate values of the parameter p (numbers near 
the curves) . In the right plot the same quantity is given as a function of p for separate values of 
mr (numbers near the curves). 



1.0 r 



0.8 



5 0.6 



A 
V 



0.2 - 



0.0^- 
0.2 




0.5 



0.4 



S 0.3 



0.2 



0.1 



0.0 



0.5 / 






0,75^ 




1 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 

P 



Figure 5: The normal vacuum stress on the boundary, (T| )b,z=o/ mD , for a minimally coupled 
scalar field in D = 4 as a function of mr (left plot) for separate values of the parameter p 
(numbers near the curves). In the right plot the same quantity is given as a function of p for 
separate values of mr (numbers near the curves). 



5 Conclusion 

In this paper we have analyzed the effects induced by a flat boundary on the VEVs of the field 
squared and the energy-momentum tensor associated with a massive scalar field in a higher- 
dimensional cosmic string spacetime. Although the analysis of the VEVs associated with quan- 
tum fields in the cosmic string spacetime taking into account the presence of boundaries, are, 
in general, developed considering cylindrical surfaces coaxial with the string, here we decided to 
adopt a different geometry by considering a flat boundary surface orthogonal to the string. The 
condition imposed on the field at the boundary is Dirichlet one. For a scalar field with Neu- 
mann boundary condition, the corresponding formulae for the VEVs of the field squared and 



16 



the energy-momentum tensor are obtained from those given above by changing the sign of the 
boundary induced parts. In the presence of the boundary, the VEVs obtained are given in terms 
of a sum of two terms: the first one due to the cosmic string itself in the absence of boundary 
and the second terms induced by the boundary. Because the analysis of the VEVs associated 
with scalar fields in a pure higher-dimensional cosmic string spacetime have been developed in 
the literature, in this paper we were more interested to investigate the contribution induced by 
the boundary. Two distinct situations were considered: first one when the parameter which 
codifies the presence of the string is the inverse of an integer number, and in the second for 
general values of this parameter. For the first case, the Green function is expressed in terms of a 
finite sum of the Macdonald functions, while for the second one, only an integral representation 
can be provided. 

For integer values of the parameter p, the boundary induced part in the VEV of the field 
squared is given by ([25]) and this part is always negative. Note that the pure topological part 
is positive and is given by the expression (I28p . For general case of the parameter p, we have 
provided the integral representation f)35[) . For points on the string axis the boundary induced 
VEV in the field squared is related to the corresponding quantity in Minkowski spacetime by 
simple formula (|36p . For points near the boundary, the leading term in the asymptotic expansion 
of the VEV does not depend on the angle deficit and it coincides with corresponding term for 
a boundary in Minkowski spacetime. At large distances from the boundary and for a massive 
field, the boundary induced part is exponentially suppressed and the VEV of the field squared 
is dominated by the pure topological part. For a massless field the boundary induced VEV of 
the field squared is given by (j38l) and at large distances it decays as |z| 2 ~-°. 

Another important local characteristic of the vacuum state is the VEV of the energy- 
momentum tensor. Similar to the field squared, the vacuum energy-momentum tensor is decom- 
posed into pure topological and boundary induced parts. For the general case of the parameter 
p, the latter is non-diagonal and is given by expression ([54"]) • For a massless field the correspond- 
ing formula is reduced to (162p . In the case p being an integer number, the general formulae take 
the forms (|64p and (|68p for massive and massless fields respectively. The corresponding VEV in 
the boundary- free cosmic string geometry is given by (|65p . For a conformally coupled massless 
scalar field the total energy density is always positive. For points near the boundary, the VEVs 
of the energy density and the diagonal stresses along the directions parallel to the boundary 
diverge as \z\~ D for a non-conformally coupled field. For a conformally coupled field the leading 
terms vanish and these VEVs behave as |z| 2_ ' D . The VEVs of the normal stress, (T| )b, and of 
the off-diagonal component, (T 3 1 )b, are finite on the boundary for points outside the string axis. 

The normal vacuum force acting on the plate is determined by the component (Tf )t, eval- 
uated on the boundary. This force is determined by expressions (|70p and (I72p for massive and 
massless fields respectively. In particular, for a massless field in 4-dimensional spacetime the cor- 
responding effective pressure is given by (j73[) and it is positive for both conformal and minimal 
couplings. Note that for a flat boundary in Minkowski spacetime the normal stress vanishes and 
the effective force in the geometry under consideration is induced by the presence of the string. 
Numerical examples presented show that the non-trivial topology due to the cosmic string en- 
hances the vacuum polarization effects for both field squared and the energy-momentum tensor 
compared to the case of a boundary in Minkowski spacetime. 

Before to finish this paper, we would like to say that the analysis developed here may be 
relevant in the investigation of vacuum polarization effects induced by a cosmic string in anti- 
de Sitter spacetime [34] . In fact for this geometry, adopting Poincare coordinate system, the 
coordinate along the string has a semi-infinity range and a flat boundary orthogonal to the string 
is present. It is also our future interest to analyze the vacuum polarization effects associated 
with quantum fields in this geometry. 
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